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We describe a simple Hamiltonian for an underdamped 
Joseplison array coupled to a single photon mode in a reso- 
nant cavity. Using a Hartree-like mean-field theory, we show 
that, for any given strength of coupling between the photon 
field and the Josephson junctions, there is a transition from 
incoherence to coherence as a function of A^, the number of 
Josephson junctions in the array. Above that value of A'^, the 
energy in the photon field is proportional to A'^^, suggestive 
of coherent emission. These features remain even when the 
junction parameters have some random variation from junc- 
tion to junction, as expected in a real array. Both of these 
features agree with recent experiments by Barbara et al. 



I. INTRODUCTION 

Researchers have long sought taxause Josephson junc- 
tion arrays to radiate coherently.BB To achieve this goal, 
a standard approach is to inject a d.c. current into an 
overdamped array. If this current is sufficiently large, 
it generates an a.c. voltage Vac across the junctions, of 
frequency loj = 2eVdc/fi, where Vdc is the time-averaged 
voltage across the junction.0 Each junction then radiates 
(typically at microwave frequencies). If the junctions are 
coherently phase-locked, the radiated power P cx N^, 
where N is the number of phase-locked junctions. This 
iV^ proportionality is a hallmark of phase coherence. But 
many difficulties inhibit phase coherence in practice. For 
example, the junctions always have a disorder-induced 
spread in critical currents, which produces a distribu- 
tion of Josephson frequencies and makes phase locking 
difficult.oEl Furthermore, in small-capacitance (and un- 
derdamped) Josephson junctionSy, quantum phase fluc- 
tuations inhibit phase locking Jj^Ej Thus, until recently, 
the most efficient coherent emission was found in two- 
dimensional arrays of overdamped Josephson junctions, 
where quantum fluctuations are minimal. 

Recently, Barbara et al. have reported a remarkable 
degree of. coherent emission in arrays of underdamped 
j unctions Ej'Ej Their arrays were placed in a microwave 
cavity, so as to couple each junction to a resonant mode 
of the cavity. If the mode has a suitable frequency and 
is coupled strongly enough to the junction, it can be ex- 
cited by a Josephson current through the junction. The 
power in this mode then feeds back into the other junc- 
tions, causing the array to phase-lock and inducing a to- 
tal power P oc iV^. For a given coupling, Barbara et al. 
found that there is a threshold number of junctions Nc 



below which no emission was observed. The coupling, 
and hence could be varied by moving the array rela- 
tive to the cavity walls. 

Barbara et al. interpreted tiiep|results by analogy with 
the Jaynes-Cummings modelE3EZl of two-level atoms in- 
teracting with a radiation fleld in a single-mode reso- 
nant cavity. In this case, each Josephson junction acts 
as a two-level atom; the coupling between the "atoms" 
is provided by the induced radiation field. A dynamical 
calculation based on a model similar to that of Jaynes 
and Cummiii gS | h as been carried out by Bonifacio and 
collaboratorst§t2l for Josephson junction arrays in a cav- 
ity. Their model does produce spontaneous emission into 
the cavity above a threshold junction number, provided 
that the Heisenberg equations are treated in a certain 
semi-classical limit appropriate to large numbers of pho- 
tons in the cavity. 

In this paper, we present a simple model for the on- 
set of phase locking and coherent emission by an under- 
damped Josephson junction array in a resonant cavity. 
We also calculate the threshold for the onset of phase co- 
herence, using a form of mean-field theory. Our model de- 
rives from more conventional models of Josephson junc- 
tion arrays, but treats the interaction with the radiation 
field quantum-mechanically. Within the mean field the- 
ory, we find that for any strength of that coupling, there 
exists a threshold number of junctions Nc in a linear array 
above which the array is coherent. Above that threshold, 
the energy in the photon field is quadratic in the number 
of junctions, as found experimentallylla The model is eas- 
ily generalized to two-dimensional arrays. Furthermore, 
as we show, the threshold condition and A^^ dependence 
of the energy in the radiation field, are both preserved 
even in the presence of the disorder which will be present 
in any realistic array. Finally, the coupling constant be- 
tween junctions and radiation field can, in principle, be 
calculated explicitly, given the geometry of the array and 
the resonant cavity. 

The remainder of this paper is organized as follows. In 
Sec.|l|, we describe our model and approximations. Our 
numerical and analytical results are presented in Sec. IB. 



Section IV presents a brief discussion and suggestions for 



future work. 



II. MODEL 
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A. Hamiltonian 

We consider a Josephson junction array containing N 
junctions arranged in series, placed in a resonant cavity, 
arranged in a geometry shown schematically in Fig. ^ It 
is assumed that there is a total time-averaged voltage $ 
across the chain of junctions; this boundary condition is 
discussed further below. The Hamiltonian for this array 
is taken as the sum of four parts: 



H = Hj + He 



phot 



int ■ 



(1) 



Here Hj = — ^Jj ^os 4>j is the Josephson coupling 

energy, where (jjj is the gauge-invariant phase difference 
across the jth junction, Ejj = hicj/q, the critical current 
of the jth Josephson junction is Icj, and q = 2\e\ is the H' = H + /i quj/Cj 

magnitude of a Cooper pair charge. He is the capacitive 
energy of the array, which we assume can be written in 



where (. . .) denotes a quantum-statistical average. We 
will impose a constant- voltage boundary condition, by re- 
quiring that $ = *i 

across the linear array should 
take on a specified value. Here $ represents the to- 
tal, time-averaged voltage across the linear array. It is 
most convenient to impose the constant-voltage bound- 
ary condition by using the method of Lagrange multi- 
pliers, adding to the Hamiltonian a term fi X) jLi ~ 

/i^^j^ qrij/Cj, where the constant /i will be determined 
later by specifying 

If we combine all these assumptions, we can finally 
write an explicit expression for H', the operator whose 
ground state we seek: 

N 



the form Hq = E^i Ecjn^j, where Ecj = q^/{2Cj), the 
capacitance of the jth junction is Cj, and nj is the dif- 
ference in the number of Cooper pairs on the two grains 
connected by the jth junction. The field energy may 
be written as iJphot ~ ?iil(a^a + 1/2), where il is the 
frequency of the cavity resonant mode (assumed to be 
the only mode supported by the cavity), and and a 
are the usual photon creation and annihilation opera- 
tors, satisfying the commutation relations [a, a^] = 1; 
[a, a] = [a''', a'''] — 0. We assume that the number 
operator rij and phase 0^ have commutation relations 
[rij, exp(±i0j;)] = ± exp(±z(/)j)^ji., which implies that Uj 
can be represented as —id/{d(j)j). 

The crucial term in the Hamiltonian for phase locking 
is the interaction term ifint- We write this in the form 
^^int = (1/c) / J- Ad'^a;, where J is the Josephson current 
density, A is the vector potential corresponding to the 
electric field of the cavity mode, c is the speed of light, 
and the integral is carried out over the cavity volume. 
Since J is comprised of the Josephson currents Icj sin ipj 
passing through the junctions, we may write this last 

N 



term as Hmt =J2-i=i ^Jj^. 



A - - 
- he 



A • ds, 



(2) 



where the integral is across the jth junction, i.e., betwewi 
the jth and (j-l-l)th superconducting grain (see Fig. |l|).E2l 
The phase factor Aj may be expressed in terms of the cre- 
ation and annihilation operators for the photon quanta 
as (in esu) Aj = iy^hc^ /{2n){a — a^)aj, where aj is a 
suitable coupling constant depending an the polarization 
and electric field of the cavity mode.cil It is convenient 
to introduce the notation hgj/^/V = Ejjaj y^hc^ /{2fl), 
where V is the cavity volume. 

Finally, we need to discuss suitable boundary condi- 
tions for this linear array. Let ^j denote the time- 
averaged voltage across the jth junction. For our as- 
sumed form of the capacitive energy, ^j = q{nj)/Cj, 



hn (a^'a + - ) + '^[^^■Jj cos(j>j + E, 



+ Hqrij/Cj H — — o)) sin0j j . 

yV J 



(3) 



B. Mean-Field Approximation 

The eigenstates of H' are many-body wave functions, 
depending on the phase variables (jij and Uj , and the pho- 
ton coordinates a and . We will estimate the ground 
state wave function and energy using a mean-field ap- 
proximation. To define this approximation, we express 
H' in the form 



H' = H 



pli; 



H 



phot 



Hu 



(4) 



where JJphaso = Y^^^ii^^Jj cos<j)j + EcjU^j + ^quj / Cj), 
and iJint = iiTi/ \/V){a — a^) X^jLi 9j sin(/)j. The mean- 
field approximation consists of writin 



AT 



Hint ~ i-i= {a - at) ^ sin 



N 



{a - a^)^gj{sm( 



N 



[a -a') 2_^gj{sm<Pj) 



(5) 



With this approximation, H' is decomposed into a sum 
of one-body terms, each of which depends only on the 
photon variables or on the phase variables of one junc- 
tion, plus a constant term. The eigenstates of H' , in 
this approximation, are of the form \E'(a, a^, {0j}) = 

'0phot(a,at) rijli V'j(<?!'i), where f/'phot and the V'j's are 
one-body wave functions. 
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That part of H' which depends on photon vari- 
ables may be written -ffp^ot = ^^phot + i{fi/y/V){a — 
a.^) "^^^i gj {sin (f>j), where (sin0j) denotes a quantum- 
mechanical expectation value with respect to 
With the definition \j = {cxp{i(j>j)) , -ffp^ot takes the form 

= nn (a^a + i) + _ at), (6) 



2/ - VV 



where Xj = {Xj — A*)/(2i) = {sin This is the Hamil- 
tonian of a displaced harmonic oscillator; its ground state 
energy eigenvalue i?phot;0 is readily found by completing 
the square to obtain 



(7) 



where fet = at_|_j ^^-^ gjXj / (fly/V), and we have defined 

(8) 



N 



(Note that b and fot have the same commutation relations 
as a and at, i. e., [b,b^ ~ 1.) 

The resulting ground state energy of H'^^iot 



Similarly in the ground state, since (6t) = 0, 



{a?) = -i r=T]. 



(9) 



(10) 



Also, the total energy stored in the photon field is 

Ep^ot = ma^a +-) = — + _7y2, (11) 



where we use Eq. (|0|) and the fact that {b^b) = in the 
ground state. 

The wave function ijjj ) is an eigenstate of the effec- 
tive single-particle Schrodinger equation Hj^pj — Ej^pj, 
where 



Hj = -Ejj COS0J -I- + 2nqnj) 



V 



ism0j(a — a'), 



(12) 



and (a — at) = 2iri/{Q\/V). Using this expression and 
completing the square, we can write 



2hgj 



nv 



rj sin <j)j — Ejj cos (pj , 



(13) 



where we have written n = —fJ./q. (Note that with our 
definition of /i, it does not have the dimension of energy.) 
Introducing the notation 



Ea;j 

where we define E\ 



Eh 



Eint;j^ 



(14) 



int;j 



tan "'^(iJintyV-E'jj), we obtain 



2hr]gj/{nV) and 



-Ea C0S( 



Ecj [{r 



The Schrodinger equation 



(15) 



(16) 



where Hj is given by expression (|T^), can be trans- 
formed into Mathieu's equation by a suitable change 
of variables. Specifically, if we use the representation 
= —id/{d<j)j), and we also make the change of vari 



ables tpj((f>j) — exp{in<l)j)uj{(l)j), then Eq. (16) takes the 
form 



{Ej 



Since (h-j and 



-Ea-j cos( 



(f>a-j)uj - Ecj 



(17) 



- 27r represent the same physical state, 
the physically significant eigenstate should satisfy 

ipj + 27r) — ipj {(f) j ) , or equivalently 

Uj{(j)j + 2tt) = cxp{~2Trin)uj{(l)j). (18) 

Thus, the solutions to Eq. ([l7| ) are Mathieu functions 
satisfying the boundary condition (|l8|). 

The total ground-state energy of the coupled system 
takes the form 



E., 



N 

E 



Ej-O + -Ephot;0 + Ed, 



(19) 



where Ej-nis the lowest eigenvalue of the Schrodinger 
equation ([l6|). Note that the i?j:o's are also functions of 
the Xj's, but only through the variable rj. Ed is a "double- 
counting correction" which compensates for the fact that 
the interaction energy is included in both -Bphot;0 and the 
Ej-o^s; it is given by the negative of the expectation value 
of the last term on the right-hand side of Eq. (f2|), i.e., 



E, = -i 



V 



N 

^'')'^9j{sin<j)j) 



2nTf 



(20) 



Hence, the total ground-state energy is 
Etotiv) 



N 



E. 



phot;0 



E, 



N 1 h 

= J2E,.o{v) + ^hn + —7j 



(21) 



The actual ground-state energy is found from this ex- 
pression by minimizing Etot with respect to the variable 
77, holding n (or n) fixed. 
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C. Approximate Minimization 

We begin by considering the case n — 0, for which 
an approximate minimization of i?tot iv) can be done an- 
alytically as follows. First, one must evaluate the en- 
ergies Ej-o, which are the ground-state eigenvalues of 
Hjipj{(f)j) :^^-Ji!jipj{(j)j). For n = 0, Ej^o has the approxi- 
mate valuecj 



Ej-Q - 

for Ea-j <C Ecj, and 
Ej-o 



2E, 



Ea:j 



(23) 



for Ea-j 3> Ecj- A function which interpolates smoothly 
between these limits is 



Er,o-Ec,~^E'c^+Ei.^^. (24) 
Substituting this expression into Eq. (pl|), wc obtain 




Setting d-Etot/dyy — 0, we obtain the condition 

V ■ 



2% V- 



E% + [2ng,/{nv)]W 



(25) 



(26) 



This equation always has the solution ry = 0. If 



2h 



N 

E 



> 1, 



(27) 



then there is also a real, nonzero solution for rj. Whenever 
this solution exists, it is a minimum in the energy, and 
the ?7 = solution is a local maximum. Thus, Eq. ( |27| ) 
represents a threshold for the onset of coherence. 

In the opposite limit, when Ea-j 3> Ecj and Xj 1, 



-Ea-j C0S( 

-i?Q,y- C0S( 



(t)a)+Ecj[{n-nf -n^] 
K)+Hl (29) 



The coherence threshold occurs in the small-coupling 
regime, \Ea-j\ ^ Ecj- The desired ground-state solu- 
tion can be obtained as a perturbation expansion about 
the solutions to the zeroth order Schrodinger equation, 
Eljipj = Ejip'j. The (unnormalized) solutions to this 
equation are 'ip'j = exp(im0j), corresponding to eigen- 

n^], with m integer. For 



(22) values E^ = Ecj[{m 



\n\ < 1/2, the ground state is to = 0. The second-order 
perturbation correction to this energy due to the pertur- 
bation H', = —Ea-j cos( 



IS 



rn— ±1 



mH',\m) 

En — E„ 



(30) 



where |m) denotes the ket corresponding to exp(zTO(/)j). 
After a little algebra, it is found that AEj = 
~El.^/[{2Ecj){l - 4n2)]. If \Ea:j\ » Ecj, then the 
ground state eigenvalue Ej-Q approaches —Ea-j as in the 



case n = 0. The generalization of the formula ( psj) to the 
case n is readily shown to be 



h ^ 



E, 




(31) 



where Ecj — Ecj{l — An^). To determine r] for a given 
value of n, and of the Ecj^s, Ejj^s, and gj's, one mini- 
mizes this energy with respect to 77, as described above. 

In practice, at any value of n, and for any given dis- 
tribution of the parameters gj, Cjj, and Ejj, one can 
easily evaluate the energy numerically, using the known 
properties of Mathieu functions, hence obtaining both 
the coherence threshold and the value of the order pa- 
rameter ry. Once 77 is known, the individual values of the 
Xj 's can obtained by numerically solving the Schrodinger 
equation (p^, using the Hamiltonian (12) for the ground- 
state eigenvalue. Finally, the constant-voltage condition 
can be imposed by choosing /i so that ^J^iq{nj)/Cj 
equals the time-averaged voltage across the array. 



N 



N 



E^J^J ^E 



9j- 



(28) 



If we define g = Y^^=i9j 1^ ^i^d = tl/^-, then we 
see that 77 rises from zero at a threshold determined 
by Eq. ( p7| ) and approaches unity when the parameters 
\Ea-j \ are sufficiently large. 

For n 7^ 0, the threshold can still be approximately 
found analytically. Since H' is periodic in n with period 
unity, one need consider only — l/2<n<l/2. In this 
regime, we write Hj as 



III. RESULTS 

Although our formalism applies equally to ordered and 
disordered arrays, we will present numerical results for 
ordered arrays only, purely for numerical convenience. 
In the ordered case, the constants gj, Ecj, and Ejj are 
independent of j. In this ordered case, we denote the 
parameters g, Eq — q^/{2C), and Ej respectively. For a 
specified value of n, we can find the ground-state eigen- 
value Ej-o numerically by solving Eq. using the well- 
known properties of the Mathieu functions. We can then 
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minimize the total energy Etot with respect to rj. In 
the ordered noted, all the A's are equal, and 

T] — NX. Furthermore, in this case, n is related to $ by 
$ — Nqn/C . Hereafter, for given values of Ec, Ej, 
and n, we define Ao as the value of A which minimizes 
the total energy i?tot ■ 

In Fig. ^, we plot Aq for this ordered array, as a func- 
tion of N , assuming n = 0. Two curves are plotted. The 
full curve shows Aq for the case Ej — 0, i.e., no direct 
Josephson coupling. The dashed curve in Fig. ^ shows 
Ao but for a finite direct Josephson coupling. In both 
cases, there is clearly a threshold array size Nc, below 
which Ao = 0. For N > Nc, we find Ao > 0. Since 
Ao — (sin0j)o (that is, the expectation value of sin^j in 
this energy- minimizing state), the Josephson array has a 
net supercurrent in this configuration. As N increases, Aq 
approaches unity, which corresponds to complete phase- 
locking. The value of N^ is larger for finite Josephson 
coupling than for zero direct coupling; thus, it appears, 
paradoxically, that the finite direct coupling actually im- 
pedes the transition to coherence. This point will be 
discussed further below. 

For Ej — and n = 0, Nc can easily be found analyt- 
ically from Eq. (p7|). The threshold is found to satisfy 

Nc^Ec/{2Ejo), (32) 

where Ejq = hg'^/{nv). This value agrees quite well 
with our numerical results (cf. Fig. ||). Note that, for 
any nonzero value of the coupling Ejq, no matter how 
small, there always exists a threshold value of N, above 
which phase coherence becomes established. 

If Ej — and n — 0, Nc can still be obtained as an 
implicit equation even in the disordered case, in terms 
of the distribution of the gj's and Ecj^s. The result is 
readily shown to be 

9n o2 

j=i -I 

For a given distribution of the parameter g^/Ecj, there 
will always exist a threshold value of N such that this 
equation is satisfied, no matter how weak the coupling 
constants gj. Thus, at least in this mean-field approxi- 
mation, the disorder has no qualitative effect on the co- 
herence transition discussed here. In particular, the crit- 
ical number Nc does not necessarily either increase or 
decrease with increasing disorder; instead, Nc depends 
on the distribution of gj , Ejj , and Ecj in the array. 

The inset to Fig. || shows the total energy in the photon 
mode, iSphot = h.^{{a^a) -\- 1/2), in the ordered case, 
plotted as a function of iV for n = 0. From Eq. (pT|), we 
find that -Ephot = + N'^XqEjq for an ordered array; 
this is the quantity plotted in the inset. As is evident 
from the plot, -Ephot varies approximately linearly with 
all the way from the coherence threshold to large 
values of iV, where Ao ^ 1. This N"^ dependence is a 
hallmark of phase coherence. 



The voltage $ across the array is determined by n (or 
equivalently fi). In Fig. ^, we plot Ao as a function of 
n for several array sizes at fixed coupling constants Ejq , 
Ec, and Ej in an ordered array. Since as already shown, 
Ao is periodic in n with a period of unity, we plot Ao (n) 
only for a single period, < n < 1. Fig._^ shows that, 
for any given N and Ejq, the calculated Ao achieves its 
maximum value when n has a half-integer value, i.e., the 
array is most easily made coherent at such values of n. 
In particular, an array whose size is slightly below the 
threshold value at integer values of n can be made to be- 
come coherent, with a nonzero Aq, when n is increased — 
that is, when a suitable voltage is applied. On the other 
hand, for values of N far above the threshold, Aq is little 
affected by a change in n. 

In Fig. 0, we show the quantity {rij) as a function of 
n, for several values of N and fixed value of the coupling 
constant ratios Ejq/Ec and Ej/Ec, for a single cycle 
(0 < n < 1) . This quantity is related to the voltage drop 
across one junction, in our model, by ^/N — q{nj)/C. 
For sufficiently large arrays, (nj) ^ n and the voltage 
drop is nearly linear in h in this mean-field approxima- 
tion. For arrays closer to the coherence threshold, (uj) is 
a highly nonlinear function of n. However, the deviation 
from linearity, {nj) — h, is, once again, a periodic func- 
tion of n with period unity. The discontinuous jumps in 
n as a function of (uj) represent regions of incoherence 
(Ao = 0), whereas the regions in which (uj) is a smooth 
function of n are regimes of phase coherence (Ao 7^ 0). 

In Fig. ^, we again plot Ao(-^) for two fixed ratios 
Ej/Ec, but this time for n = 1/2. From Fig. ||, we ex- 
pect this choice of n to maximize the tendency to phase 
coherence and thus to reduce the threshold array size for 
the onset of phase coherence. Indeed, in the absence of 
direct Josephson coupling, this threshold is reduced to 
below unity (that is, Xq remains nonzero, even at = 1, 
for our choice of Ejq). In fact, for this value of n, only 
an infinitesimal coupling to the resonant mode is required 
to induce phase coherence in this model. Once again (cf. 
Fig. ^) , the addition of a finite direct Josephson coupling 
actually increases the threshold number for phase coher- 
ence at n = 1/2 as it does at n = 0. 

Although we have not carried out a similar series of 
calculations for a disordered, our analytical results show 
that the essential features found in the ordered case will 
be preserved also in a disordered array. Most impor- 
tantly, there remains a critical junction number for phase 
coherence in a disordered array, just as there does in the 
ordered case. The most important difference between the 
two cases is that the individual Aj's will be functions of 
j in the disordered case. 

IV. DISCUSSION 

Although the present work is only a mean-field approx- 
imation, we expect that it will be quite accurate for large 
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N. The reason is that, in this model, the one photonic 
degree of freedom is coupled to every phase difference, 
and thus experiences an environment which is very close 
to the mean, whatever the state of the individual junc- 
tions. Such small fluctuations are necessary in order for 
a mean-field approach to be accurate. In fact, a simi- 
lar approach has proven very successful in work on novel 
Josephson arrays in which each wire is coupled to. aJarge 
number of other wires via Josephson tunneling.cj'c3 

It may appear surprising that a finite direct Joseph- 
son coupling actually increases the threshold array size 
for coherence. But in fact this behavior is reasonable. If 
there is no direct couphng {Ej = 0), the phase difference 
across each junction evolves independently, except for the 
global coupling to the resonant photon mode. When the 
array exceeds its critical size, this coupling produces co- 
herence. If the same array now has a finite Ej, there are 
two coupling terms. But these are not simply additive, 
but in fact are 7r/2 out of phase: the direct coupling fa- 
vors (f>j — 0, while the photonic one favors <j)j = Tr/2. For 
a large enough array, the coupling to the photon field still 
predominates and produces global phase coherence, but 
this occurs at a higher threshold, at least in our model, 
than in the absence of direct coupling. 

A striking feature of our results is the very low coher- 
ence threshold (A'^ = 1) when n = 1/2. In fact, for any 
N and for Ej = 0, only an infinitesimal coupling to the 
cavity mode would be required to induce phase coherence 
at n = 1/2. The reason for this low threshold is that, in 
the absence of coupling, junction states with {rij) = 
and (rij) — 1 are degenerate. Any coupling is therefore 
sufficient to break the degeneracy and produce phase co- 
herence. A related effect has been noted previously in 
studies of more conventional Josepihaon junction arrays 
in the presence of an offset voltage.t3 

Finally, we comment on what is not included in the 
present work. This paper really considers only the mini- 
mum energy state of the coupled photon/junction array 
system under the assumption that a particular voltage is 
applied across the array. It would be of equal or greater 
interest to consider the dynamical response of such an 
array. Specifically, it would be valuable to develop and 
solve a set of coupled dynamical equations which incorpo- 
rate both the junction and the photonic degrees of free- 
dom. Such a set of equations has already been proposed 
by Bonifacio under a particular set of simplifying assump- 
tions. A more accurate set of equations is needed, which 
would include not only a driving current, but also the 
damping arising from both resistive losses in the junc- 
tions and losses due to the finite Q of the cavity. In the 
absence of damping, such equations can be written down 
from the Heisenberg equations of motion. The inclusion 
of damping may be more difficult. We hope to discuss 
some of these effects in a future publication. 
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FIG. 5. Same as Fig. |, but for n = 1/2. 



FIG. 1. Schematic of the geometry used in our calculations, 
consisting of an underdamped array of Josephson junctions 
coupled to a resonant cavity, and subjected to an applied 
voltage $. The array consists of + 1 grains, represented 
by the dots, coupled together by Josephson junctions, repre- 
sented by the crosses. Nearest-neighbor grains j and j + 1 
are connected by a Josephson junction, and the cavity is as- 
sumed to support a single resonant photonic mode. For the 
specific calculations carried out in this paper, we assume a 
one- dimensional array as shown, and a specific form for the 
capacitive energy, as discussed in the text. 



FIG. 2. Coherence order parameter Ao which minimizes 
£^tot(A) for a one-dimensional array, plotted as a function 
of the number of junctions A'^, for two values of the di- 
rect Josephson coupling energy Ej. Other parameters are 
hg/VV = O.SEc, hn/2 = 2.6Ec, and n = 0. The coupling 
parameter Ejo = hg^/(QV) is given by Ejo ~ 0.017i5c- In- 
set: total energy in the photon field, -Ephot, plotted as a func- 
tion of A^, for the same parameters and the same two values 
of Ej. 



FIG. 3. Energy-minimizing value of the coherence order 
parameter Aq, as a function of the parameter n = n/q, 
for several values of the array size N. Other parame- 
ters are hg/VV = O.SEc, hQ/2 = 2.6Ec, Ej = 0, and 
Ejo X 0.017 Ec. 



FIG. 4. The parameter (rij) = $C/{qN), where <1>/A is 
the voltage drop across one junction, plotted as a function of 
the parameter n — ^/q, for several values of the array size 
N. Other parameters are hg/VV = 0.3£'c, fi^/2 = 2.6Ec, 
Ej = 0, and Ejo « 0.017-Ec. 
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Figure 1. 

Model for a Josephson junction array coupled to a resonant cavity. 
J. Kent Harbaugh and D. Stroud. 
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Model for a Josephson junction array coupled to a resonant cavity. 
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Figure 4. 

Model for a Josephson junction array coupled to a resonant cavity. 
J. Kent Harbaugh and D. Stroud. 
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